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We study a class of homogeneous finite-dimensional Ising models which were recently shown to 
exhibit glassy properties. Monte Carlo simulations of a particular three-dimensional model in this 
class show that the glassy phase obtained under slow cooling is dominated by large scale excitations 
whose energy Ei scales with their size Z as i?; ~ with Q 1.33(5). Simulations suggest that 
in another model of this class, namely the four-spin model, energy is concentrated mainly in linear 
defects making also in this case domain walls tensionless. Two-dimensinal variants of these models 
are trivial and energy of excitations scales with the exponent — 1.05(5). 



Recently the problem of the structure of the glassy phase in spin glasses has attracted considerable attention. The 
main thrust of research has been to establish whether the low-temperature phase is described by the so-called " droplet 
model" or by the replica symmetry breaking (RSB) theory of Parisi g]. An important difference between these 
theories concerns the energy Ei of large-scale excitations which should scale with their size I as Ei ^ where > for 
the droplet model but with = for the RSB theory. The most interesting situation arises in the three-dimensional 
case, where it seems that a combination of these two approaches is needed to describe the glassy phase correctly [||. 

Of course, the problem of the structure of the glassy phase is not restricted to spin glasses. Disordered, out of 
equilibrium slow dynamics structures appear in superconducting compounds, polymers and granular matter. However, 
these systems are very complex and modelling them seems to be more difficult than understanding spin glasses. 
Conventional glasses are also very complex Nevertheless, recently relatively simple models have been proposed 
which exhibit an encouraging number of glassy properties ■ AH these models are spin models which do not contain 
quenched disorder as in the case of spin glasses and glassiness is dynamically generated. An absence of quenched 
disorder has important implications. Firstly, their ground state and sometimes even the structure of excitations are 
known exactly. Let us emphasize that for the three-dimensional spin glasses the problem of finding the ground state is 
extremely difficult (NP-complete) and is one of the main difhculty in numerical approaches to spin glasses. Secondly, 
for models without quenched disorder there is no need to average over different realizations of this disorder, which is 
yet another computational advantage of such models. 

The objective of the present paper is to examine the nature of glassy phase in certain nondisordered Ising models. 
Simple heuristic arguments show that these models might have large scale excitations of energy which scale with their 
size as i.e., slower than their surface (~ The question is whether such states appear in, or maybe even 

dominate, the glassy phase. To examine this problem we have performed Monte Carlo simulations of the models and 
our results suggest that the glassy phase is dominated by excitations whose energy increases faster than /^^^^ but 
slower than l'^^^. It is also likely that the exponent 8 which describes the size dependence of energy of excitations 
takes some nontrivial values for these models. For the models of conventional glasses considered here the problem of 
energetics of large scale excitations is computationally much more tractable than for spin glass models. Hopefully, 
results obtained for these models will provide some insight into other glassy systems too. In addition, our results can 
be used to verify some earlier claims concerning the nature of the glassy transition in some of these models . 

The class of models which we examine in this paper is defined by the following Hamiltonian 

H=-2kY, S^S, + ^ ^^^1 + ^^-^ E S^SJSkSl, (1) 

<i,j> «i.J>> li,j,k,l] 

where summations in (|l|) are over nearest neighbours, next-nearest neighbours, and elementary plaquettes respectively. 
Model (|^) has the interesting property that the energy of cerain excitations of size I is proportional to Z'^"^, whereas 
typically the energy of an excitation for a standard nearest neighbour Ising model is proportional to its surface 
(~ This property has been used to construct a class of random surface theories based on model (|l|) pl] |. 

Recently, it was shown that for d = 3 model ([|) has slow dynamics at low temperature (lO| . When the high temper- 
ature sample is quenched to low temperature the excess energy SE = E(t) — Egg, where E^q is the equilibrium energy, 
decays with time t much more slowly than which is a typical decay rate for Ising models with nonconservative 

dynamics p2| . It turns out that the fc = (pure four-spin interaction) case is of particular interest. This is because in 
this case the model has also some other properties typical of conventional glasses such as strong metastability ||^ and 
small cooling-rate effects p^ ]. Moreover, certain time dependent correlation functions, such as those describing aging, 
also behave similarly to real glassy systems p3. Although a slow decay of SE is an indication of slow dynamics, it 
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would be desirable to relate this decay with the increase of a characteristic length scale I. (As we will see such a 
relation will give some information about the energetis of excitations of the glassy phase.) 

For ordinary Ising models simple arguments, based on the fact that the energy of excitation of the size / scale as 
its surface (£"/ ^ ^'^~^), can be used to obtain the relation 

SE^j. (2) 

However, since for model energy of excitations might increase more slowly than their surface area, the relation (^) is 
no longer obvious. Assuming that in the glassy phase the dominant excitations are indeed these low-energy excitations 
(with El - /'^"^), the following relation should hold 

<5i?~l. (3) 

Let us notice that the assumption that the glassy phase is composed of low-energy excitations implies that at the 
glassy transition domain walls lose their surface tension. Such an identification might be of more general validity and 
could be used as a new criterion to locate the glassy transition. 

Is it possible to verify which of the relations (||) and (||) are true in our model? First, let us notice that (||) and (|^) 
are two extremal cases corresponding to the largest and the smallest excitation energy per surface area, respectively. 
It is thus possible that neither of them is true and in the glassy phase an intermediate relation holds. To consider a 
more general situation let us assume that energy of excitations scale as . In a lattice of the linear size L the number 
of excitations of size I scales as {L/lY and the total excess energy scales as [L/lYl® . Thus, the excess energy per 
spin 6E scales as 

5E ~ l^-"^. (4) 

To determine Q we need a second, independent measurement of the characteristic length /. It is already known 
that I can be also obtained from the fluctuations of the order parameter using the relation [|5| 

X = \/L'' <llL\Y.S,f>^l\ (5) 

i 

where the magnetization is taken as a corresponding order parameter. Assuming that relation (^ determines the 
same characteristic length as (^ (up to the order of magnitude) , we use the following procedure: We continuously cool 
the high-temperature sample down to zero temperature. In this process temperature changes linearly with simulation 
time 

T{t) = To - Tt, (6) 

where r is the cooling rate. Then, for the zero-temperature configuration we calculate 5E and x- Previous Monte 
Carlo simulations suggest that for model (|I]) and fc = 2 the zero-temperature energy excess 6E decreases with r as 

6Er^r-'\ (7) 

where xi = 0.50(5) [ p^ . Similarly, we expect that x also increases as a power of r 

X-r--'-- (8) 

Inverting (||) and using (^) we obtain r ~ l^'^/'-^^ and from (0) we have 5E ~ i-dxi/x2_ Yin&Wy, comparing the last 
relation with (^) we obtain 

Q = d{l^xi/x2). (9) 

To estimate we made Monte Carlo simulations of model (|l|) for k = 2 and using a sequential Metropolis algo- 
rithm Simulations were made for the system size up to i = 70 and we have checked that this is sufficient to 
obtain basically size independent results. For each cooling rate r we made around 100 independent runs which were 
used to calculate 6E and < l/L'^Yl >• The starting temperature was T = 2.8, which for fc = 2 is above 
the critical temperature, which in this case is ^ 2.35 [ p^ . 

The results of our simulations are shown in Fig. |l|. The relatively good linearity of our data confirm power-law 
behaviour (|^) and (||). From this data we estimate xi = 0.50(5), X2 = 0.90(5) and using (||) we obtain 6 = 1.33(5). 
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Such an estimate of Q shows that neither (g), which corresponds to Q = 2 nor (Q), which corresponds to © = 1 are 
correct. Instead, we have an intermediate possibihty with a noninteger value of Q. Let us notice that since 6 < 2 the 
surface tension of domain walls vanishes. 

For comparison, in Fig. ^ we also present results of our simulations for the two-dimensional (square lattice) version 
of model (|l|) with fc = 2 (in this case Tc — 0). Simulations were made for system sizes up to L = 1000. From this data 
we estimate Xi = 0.46(5), X2 — 0.95(5) and thus = 1.05(5). It is likely that in this case 8 = 1, which would indicate 
a trivial nature of domain walls with energy proportional to their perimeter, the typical 2D Ising model behaviour. 

Although for k — 2 model (^ has slow low-temperature dynamics, it does not display a genuine glassy transition. 
As we already noted, to model glassy transitions one should really study the case of fc = 0. However, in this case the 
above method encounters some difficulties since domains of random quench are not only of ferromagnetic type as in 
the case of fc = 2, but also antifcrromagnetic and even of some mixed types (see |l^ for some discussion). For fc = 
equation (^ cannot be used and thus, the exponent cannot be determined using the above method. 

To get some insight into the fc = case we instead looked at the distribution of unsatisfied plaquettes in the 
glassy phase (i.e., plaquettes contributing energy above the ground state). The random high-temperature sample was 
slowly cooled down to zero temperature. Then for the final configuration we located unsatisfied plaquettes and their 
spatial distribution is shown in Fig. |^. For comparison we also present similar calculations for the k — 2 case. One can 
see that in both cases energy is concentrated in linear segments. For fc = 2 this is in agreement with our estimation 
< 2 as for Q = 2 energy would be localized on two-dimensional surfaces. Although for fc = we cannot estimate 
O, the linear stuctures in Fig. || strongly suggests that in this case also 9 < 2 and the glassy phase is composed of 
tensionless domain walls. 

Of course, the glassy phase obtained by the finite-rate cooling contains some regions other than linear segments 
where energy is concentrated. However, we expect that such spots will diminish for decreasing cooling rate r. To 
confirm our expectations we measured the ratio p of unsatisfied plaquettes which belong to linear segments compared 
to the total number of unsatisfied plaquettes |l8|. The results, presented in Fig. ^, show that p indeed increases for 
decreasing r. It is also possible that in the limit r — > the fraction p ^ 1. Let us notice that a glassy state obtained 
in such a limit constitutes an ideal glass |^9| and the present results might shed some light on this, so far hypothetical, 
state of matter. In particular they suggest that in the ideal glass the slow cooling removes energy-rich spots and 
leaves only low-energy excitations. 

In conclusion, we studied the glassy phase of gonihedric model. Our result show that energy of excitations in this 
phase scales as with Q < d — 1. It confirms earlier expectations that domain walls in this model are tensionless. 
Since the fc = case seems to have a number of properties typical to ordinary glasses, it would be desirable to check 
whether this result has also some analogy in real systems. 
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FIG. 1. The excess energy 5E ( +(d = 3) and A{d — 2)), and x ^ {*{d ~ 3) and x(d = 2)) as a function of the cooling 
rate r. 
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k=0, r=0.00001 



k=2, r=0.0005 




FIG. 2. A distribution of unsatisfied plaquettes in the zero-temperature glassy phase. Simulations were made for the system 
size L — 100. (Only a portion of the system is shown.) Let us notice that although the cooling for the k — 2 case is faster it 
seems to create larger domains. It was already suggested that for the case k = the model should order much more slowly 
than for fc = 2 0. 
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FIG. 3. The fraction p of unsatisfied plaquettes which are located in linear segments as a function of the cooling rate r for 
fc = (+) and fc = 2 (x). 
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